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The collapse of flexible polyelectrolytes in a solution of multivalent counterions is studied by 
means of a two state model. The states correspond to rod-like and spherically collapsed conforma- 
tions respectively. We focus on the very dilute monomer concentration regime where the collapse 
transition is found to occur when the charge of the multivalent salt is comparable (but smaller) to 
that of the monomers. The main contribution to the free energy of the collapsed conformation is 
linear in the number of monomers A'^, since the internal state of the collapsed polymer approaches 
that of an amorphous ionic solid. The free energy of the rod-like state grows as A*' In A'', due to the 
electrostatic energy associated with that shape. We show that practically all multivalent counterions 
added to the system are condensed into the polymer chain, even before the collapse. 



I. INTRODUCTION 

Nucleic acids and many proteins are charged and ex- 
hibit self-assembly. For example, DNA with an extended 
length of the order of centimeters is highly compacted in 
cells and viruses. The compactness of strongly charged 
chains in low ionic strength solutions has been a sub- 
ject of recent great interest Strongly charged 
chains, though in principle in extended conformations to 
decrease their electrostatic energy, precipitate into highly 
compact structures when small amounts of high valence 
counterions are added to low ionic strength solutions 
P,|6|p^. Flexible polymers, including single stranded 
DNA and many synthetic polymers such as Polyester 
Sulphonated, precipitate into dense amorphous struc- 
tures Semi-flexible polymers, such as long double 
stranded DNA, and rigid rods, such as short DNA frag- 
ments, aggregate into highly compacted ordered struc- 
tures Multivalent induced DNA compaction is a 
promising mechanism to "pack" DNA, a critical problem 
in gene therapy given that long extended DNA cannot be 
deliver into cells. Furthermore, the precipitation is corre- 
lated with highly accelerated rates of DNA renaturation 
and cyclization ||l^ , [l5|| , and other important biotechno- 
logical processes. Aside from the DNA compaction and 
its importance in biology and biotechnology, the physics 
of self-assembly of charged molecules remains an inter- 
esting open problem. 

The self-assembly of polyelectrolytes is strongly re- 
lated to the "condensation" of the counterions along 
the chains. A fraction of the counterions will be con- 
densed, strongly attached to the backbone of the poly- 
electrolytes to decrease their electrostatic energy. In flex- 
ible chains the interaction between the polyelectrolytes 
and the counterions is non-trivial. The long range elec- 
trostatic interaction that leads to ion condensation is 
strongly dependent on the chain conformation; more ions 
are condensed when the chain is compacted than when it 
is extended. The conformation of the chain, reciprocally, 
depends on the number of condensed counterions, and 
their distribution along the chain. The correlations be- 



tween monomers and counterions are such that the higher 
the number of condensed counterions, the smaller the 
size of the chain This mechanism induces a transi- 
tion between two possible polyelectrolyte conformations: 
a highly compact shape with nearly zero effective charge 
(it is neutralized by the condensed counterions) and a 
rod like shape with a smaller number of condensed coun- 
terions. 

The above counterion mediated precipitation mecha- 
nism is quite general. However, it is non trivial to com- 
pute the free energy of the collapsed chains. In this paper 
we use simple physical arguments and known results from 
related charged systems to construct the free energy of 
the collapsed state of flexible polyelectrolytes. We aim to 
describe the multivalent induced precipitation in highly 
dilute polyelectrolyte solutions at low ionic strength so- 
lutions of monovalent ions. It is well documented exper- 
imentally 1^,^ that in this regime the polyelectrolytes 
precipitate when the concentration of high valence ions 
is proportional to the concentration of charged monomers 
and that the transition is quite universal in this regime as 
it is independent of the charge density of the chains and 
on the specific nature of the condensing agents. In this 
paper we recover this result with a first order transition 
model from rod-like to collapsed states described below. 
We find that at the precipitation point all the multivalent 
ions replace the condensed monovalent ions, as expected, 
and that the effective charge of the finite size precipitate 
is nearly zero. 

II. MODELING CONSIDERATIONS 

In our approach we assume that there are two well 
defined types of states allowed for the single chain: col- 
lapsed and extended. We calculate free energies for each 
of these types of states and study the conditions under 
which one of the states is preferred. Ideally, we would 
like to present a single free energy functional valid for the 
whole range of conformations of the chain. It is, however, 
very hard to construct such a functional. On the other 
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hand, as we aim to show below, there are very distinct 
characteristics in the two phases so that treating them as 
different objects is not only possible, but provides a clear 
view of the reasons behind the transition from extended 
to collapsed conformations. 

It has long been understood how the extended states 
arise. The bare charge of a polyelectrolyte is always re- 
duced by the condensed counterions, but it is neverthe- 
less sufficient to create repulsions between the segments 
of the chains that overcome the entropic elasticity of the 
chain, leading to an almost fully stretched conformation. 
The collapsed state is, comparatively, much harder to un- 
derstand. The repulsive interaction between the uncom- 
pensated charges is still present, and therefore, a suffi- 
ciently strong cohesive energy is necessary to maintain 
the collapsed state. This cohesive energy has been mod- 
eled in previous works by means of virial approximations 
|0,D and other averages over positional fluctuations of 
the charged particles in the system We pro- 

pose here a different approach to the calculation of the 
cohesive energy for the collapsed state . 

In several treatments of the problem of attractions 
between rod-like polyelectrolytes the attractive energies 
have been identified as arising from strong correlations 
between the condensed counterions. These correlations 
can be so strong as to induce the formation of crystalline 
order of the condensed counterions. In the same way, 
we argue that also in the case of flexible polyelectrolytes, 
purely electrostatic interactions can lead to the collapse 
of the chains. A second degree of difficulty arises from the 
fact that even after identifying the origin of the attrac- 
tive interactions, the calculation of the precise value of 
the energy of interaction can become a technically diffi- 
cult problem. In this paper we will adopt the view that in 
the collapsed state, the monomers and condensed coun- 
terions of a sufficiently large polymer approach a homo- 
geneous bulk state inside the collapsed region. There, 
we can then estimate the interaction energies by consid- 
ering a similarly composed amorphous ionic solid. One 
lesson learned from the study of association of rigid poly- 
electrolytes is that it is possible to estimate interaction 
energies by simply using low temperature limits 
That is, within the collapsed region, the entropic consid- 
erations produce only small corrections. 

We can indeed consider a thought-experiment in which 
we cut the polymer bonds so as to have a mixture of 
single-charge monomers and counterions in solution. In 
this case we expect to observe the typical behavior of a 
salt. When the amount of solvent is large, all particles 
are dissociated, and as the solvent amount decreases, we 
can observe the formation of salt crystals. This process 
occurs at common temperatures and concentrations. Re- 
joining the monomers to form again the polymer chain 
modifies the previous picture in two ways. First, the dis- 
sociation in the presence of large amounts of solvent can 
never be complete as the counterions can be freed but not 



the monomeric charges. Second, in the collapsed state, 
at low solvent concentrations, the crystallization is hin- 
dered by the connectivity of the chain, which provides 
an entropic barrier to full crystallization and imposes a 
clear, well defined size for the minimum size of the crystal 
grain (the collapsed radius of gyration). 

On the basis of these considerations we obtain a model 
free energy for both states, and calculate the amount 
of added multivalent salt required to collapse the chain. 
While the chain is expected to collapse even without the 
addition of multivalent salt, such collapse would occur 
only at very low temperatures, and much of the interest 
on this process is related to applications at room temper- 
atures. The reason why the multivalent counterions are 
so effective at collapsing the chain also becomes apparent 
in the elucidation of the cohesive energy of the collapsed 
chain. Since this energy is dominated by the electrostatic 
interaction between particles, the replacement of mono- 
valent by multivalent charges has the effect of highly de- 
creasing the self-energy part of the counterions. This is 
further helped by the fact that it is entropically much 
more favorable to free multivalent instead of monovalent 
counterions. As we will show below, this replacement 
of charges is almost complete even in the non-collapsed 
state. 

III. FREE ENERGY CALCULATIONS 

We describe now the basic parameters of the system. 
We consider monodisperse polyelectrolyte chains with 
monomer number N that we take to coincide with the 
number of monovalent, positive charges in the chain. We 
assume that the distance along the chain between charges 
is a constant b that is comparable with all other small 
distances in the system such as the thickness (diameter) 
of the chain and the non-electrostatic persistence length. 
The monomer concentration is (j), and the solution con- 
tains several other charged components: a monovalent 
salt with concentration s, and a multivalent salt with 
concentration m, and valence z. We assume that the 
monovalent counterions naturally associated to the chain 
are identical with the dissociated negative components of 
the monovalent salt. The concentration of these negative 
ions is then s + (p. The size of all the small particles is also 
taken to be order b. The dimensionless Bjerrum number 
is the ratio Bj = / AnebkBT , with e the electron charge, 
e the permitivity of water, ks the Boltzmann constant, 
and T the temperature. 

A. Energy of free charges 

Whether the chain is in a collapsed or in an extended 
state, the charges that are not condensed to the chain 
give rise to similar contribution to the free energy mostly 
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arising form their translational entropy, and we will con- 
sider them all at once. A fraction / of the total charge 
in every chain is compensated by condensed counterions, 
that remain in close proximity to the chain. This com- 
pensating charge is made up from monovalent and multi- 
valent counterions. The respective charge fractions that 
they compensate are /s, and /„, so that / = fs + fm- 
Per chain, we have fsN, and f^N/z, condensed nega- 
tive monovalent and multivalent ions per chain. The free 
positive ions per chain of the two different species are 
simply sN/4>, and mN/(j)^ since none of the positive ions 
condense to the chain. Finally, we have (1 — fs)N + s/cj) 
free negative monovalent ions and mN / (j) — fmN / z free 
negative multivalent ions. 

For sufhciently dilute concentrations of polymers the 
electrostatic contribution of the free charges can be 
shown to be very small. In a Debye-Huckel approxima- 
tion, the electrostatic energy (in a per chain basis) is 
proportional to (j)^/'^. Instead, the main entropic contri- 
bution is proportional to ln0. The volume occupied by 
the chain is considered negligible and we obtain for the 
translational entropy of the free ions, per monomer, in 
units of ksT: 

F,^l Inb^s + - frn/z) Inborn - /™0/z) 

+ ^ Inborn fs) Inb^is + (1 - /«)</)). (1) 

In this expression we have estimated the entropy by 
counting the number of sites available to the ions in a 
lattice with sites of size b that are always occupied by 
ions or solvent molecules. Since the system is dilute we 
have not included the terms of the form ^(1 — s)b^ etc. 



B. The Extended State 

To estimate the electrostatic free energy for the ex- 
tended state, we assume a simple linear geometry for the 
chain. We assume that the chain is strongly stretched so 
that we take the linear size to be simply L = Nb. The 
elastic entropy of the chain plays only a small role, and 
reduces the effective size only by a small fraction. Assum- 
ing a uniform distribution of condensed charges along the 
chain, the effective charge per monomer is 1 — /s — /,„. 
The main contribution to the electrostatic free energy 
arises from the interaction of the charges with the loga- 
rithmic potential created by the effective charge, 

Fe,=Bjil-fs~f^flnN. (2) 

For the purpose of comparison with the collapsed state, 
we note that this energy is quadratic in the fraction con- 
densed, and dependent on the size of the system. 



C. The collapsed state 

We begin with the assumption that the collapsed state 
is very dense. Taking into account the space needed 
to accommodate the condensed counterions, we expect 
the radius of the collapsed polymer to be roughly R = 
{2Nb^y^^. To calculate the electrostatic free energy 
of the system we need to determine the distribution of 
charges within the collapsed polymer. To model this dis- 
tribution we borrow two arguments from classical elec- 
trodynamics and solid state physics. 

First, we consider the distribution of uncompensated 
charges. In the classical (i.e. non-statistical) analog 
of this problem, whenever there are free mobile charges 
within a solid of fixed shape, the excess charge resides 
in the surface of the object. Thus we can argue that 
the counterions in the collapsed polymer occupy places 
"deep" , in the collapsed region to cancel the monomeric 
charge there, while the few uncompensated monomeric 
charges will appear at the surface of the ball. The ther- 
mal fluctuations of the system decrease the sharpness of 
this result, and we expect a broadening of the regions in 
which there can be detected a non-zero average charge. 
We separate the collapsed ball into these two regions, 
the charged surface and the neutral core. Assuming a 
spherical geometry for the collapsed region, the interac- 
tion of the charges at the surface is readily calculated by 
an application of Gauss' law, and is simply 

Fe, = B,(l-/.-/e)^^^^. (3) 

As we shall see later, the fraction condensed for the col- 
lapsed state is very close to unity, and the energetic con- 
tribution to this charged shell is very small. 

Let us consider now the self-energy of the core, which 
is, in the average, electrically neutral. When the poly- 
mer is collapsed, all charges are very close to each other, 
and their electrostatic interactions are very strong. For 
example, for the case of the aggregation of rod-like poly- 
electrolytes, several model calculations have shown that 
there are very strong correlations between the condensed 
counterions leading to the formation of Wigner crystals 
along the rods. In the present case, however, we cannot 
argue for such crystallization because of the constraints 
that the connectivity of the chain imposes. Instead, we 
propose to view the core as an amorphous ionic solid. Fi- 
nally, while we do not expect to find a crystalline state, 
and only observe very strong correlations, the evalua- 
tion of the free energy can be carried out in a very simi- 
lar manner to ionic crystals. This amorphous ionic solid 
state has also been proposed recently by Lee and Thiru- 
malai that they call Wigner Glass |l|. 

In the core there is a very strong energetic penalty for 
conformations in which the charge of a particle is not 
quickly compensated by its neighbors. This is clearly 
the case as well in fully crystalline structures, but this 
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requirement is essentially all that is needed for a sim- 
ple estimation of the electrostatic bulk energy of crystal, 
and therefore is also applicable to our case. The energy 
per atom in an ionic crystal is usually presented in the 
following form: 



Eeb — 



(4) 



where is the lattice constant and M is the Madelung 
constant The Madelung constant depends on the 

specific lattice structure, but its value is always of or- 
der 1. The sum of the interaction of one atom with all 
other atoms in the crystal is arranged in a convergent 
series of alternating terms, summing the interaction with 
a layer of oppositely charged atoms, (nearest neighbors) 
then with the layer of next-nearest neighbors that are 
equally charged, etc. Clearly, the final result is dom- 
inated by the interaction with the neighbors that com- 
pensate the charge of a given particle, and the rest simply 
modify the prefactor. In the following we sketch a sim- 
ple estimate for an effective Madelung constant for the 
collapsed polymer bulk. 

We have assumed that all small distances in the sys- 
tem are of order b. We cluster one counterion with all the 
monomers that are required to compensate its charge, 1 
for monovalent, and z for multivalent ions. If all the par- 
ticles in the cluster are at about the same distance from 
each other, the energy contribution of the cluster is 



kT -'2^ ' 



(5) 



This is the first term in the series for the effective 
Madelung constant, and we will truncate the series at this 
point. There are fsN clusters with monovalent counte- 
rions, and fmN/z with multivalent counterions, so that 
our estimate for the free energy of the bulk of the col- 
lapsed state is 
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IV. PHASE TRANSITION CONDITIONS. 

We can now use the free energy expressions postulated 
for the collapsed and extended states to calculate the 
phase diagram of the system which we shall present in 
(j) — m coordinates. It is useful, however, to carry out 
further simplifications in the model that lead to some 
qualitative and quantitative results. 

A. Replacement by multivalent counterions 

The general expressions for the free energy in the two 
states, collapsed: 



and extended: 



Fr = Fpi 



F ^ F 



F, 



F,, 



(7) 
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still depend on the fractions condensed fs,fm- These 
fractions should be obtained before computing the final 
free energy. To obtain their values, the free energy is 
minimized with respect to these variables, which leads to 
the equations: 



0, 



2Bj{l - fm - fs) IniV + {I - fs)(b) 

2Bj{\ - ,U - fs) lnN+ (1/z) \nb^{m - f„,4>/z) = 0. (9) 

for the extended state, while for the collapsed state one 
has: 



1^(1 - fs - fm) + Bj + \nb^{s + (1 - fs)4>) = 0, 
B]N 

iVl/3 



%il-fs-f.) + Bf-±l 



iln63(m-^)=0. (10) 



The Eqs. (^ imply a particular relation between the 
fraction condensed of the two species, regardless of the 
strength of the electrostatic term (measured by Bj), 
given by 



(53(s + (l-/,)0)- =53(m-/„</./z). 



(11) 



This relation is essentially a result of the fact the same 
energy binding z monovalent counterions is what is re- 
quired to bind one multivalent ion, while the entropy 
gain of freeing z ions is also z times bigger. For dilute 
concentrations of monomers and ions, the expressions in 
Eq. (jl^) are all small quantities and even for z = 2 the 
exponentiation of the left hand side of the equation leads 
to an extremely small value for the difference m — fm4>l 
Thus we conclude that for all practical purposes all mul- 
tivalent counterions are condensed. At higher concentra- 
tion of these counterions, there can be a near saturation 
since the polymer cannot absorb all of them, and it can- 
not be emptied of monovalent ions (since by Eq. (Ill]) fs 
cannot be zero.) 

For the second pair of equations Eq. |l^ we note that 
the first terms of the equations, which arise from the 
electrostatic repulsion of uncompensated charges, present 
very large prefactors to the effective charge 1 — /s — fm- 
Clearly, the effective charge must cancel this contribu- 
tion, and we conclude that, essentially, in the collapsed 
state all charges are compensated. An analysis simi- 
lar to the case of the extended conformation implies 
that the charge compensation must occur first with as 
many multivalent ions as possible, and then finalize with 
monovalents so that, roughly, we have fm = zm/4>, and 

fs~^ fm- 



4 



B. Low Density Limit 



C. Dense limit 



Let us now consider the case in which the density of 
monomers is very low. More precisely we define this 
regime by the condition that if all chains are extended, 
they do not overlap. Since the extended size of a chain 
is Nb, the definition is equivalent to </> < N/{Nb)^ = 
1/N%3. 

To obtain a simplified formula for the amount of mul- 
tivalent counterions needed to produce the transition, we 
need to carry out some simplifications. As the solution 
becomes more and more dilute, there is a great entropic 
energy to be gained by freeing the monovalent condensed 
ions. At the same time, the relation obtained between the 
amount of condensed ions of the two species still implies 
that a large number of multivalent counterions are con- 
densed. Thus, a first simplification consist on setting, 
for the extended state fs = 0, and fm = zm/4>. Thus, 
neglecting some other small terms, the energy for the 
extended state reads: 

)' In iV-(l + -) In 6^(0 + 5). (12) 
<P 9 

In the collapsed state, we have already discussed that at 
all concentrations we have almost exactly /„ = mz/0, 
and /s — 1 — /.,„, which leads to 



-Bj{l 



z{z — l)m zm 
20 ' ~ ~ 



In 6^ 



^ + s). (13) 



In the evaluation of the entropic term we have used the 
fact that the density of required counterions is of the 
order of the monomer density, so that hib^fs{<j) + s) « 
In b^{(j) + s). 

Equating the previous two approximations for the free 
energies in the two states leads to a quadratic equation in 
zm for the density of multivalent counterions required to 
create the transition. However, from the analysis of these 
equations, a simpler conclusion can be drawn. By looking 
to just the entropic type terms, of the form ln6'^(0 -f s), 
which are the largest contributions to these expressions, 
it is clear that a transition is only possible when the 
charge of the multivalent particles approaches that of the 



monomers: 



(14) 



This is the resulting limiting condition for very dilute so- 
lutions. We should remark that while this result might 
appear intuitive, it relies on the specific structure of the 
extended and collapsed state. Namely, in both states, 
all multivalent ions are trapped, in the collapses state, 
the only free monovalent counterions are those that have 
been replaced by multivalent ones, and in the extended 
state, all monovalent counterions are free. 



In the opposite case, for higher concentrations, the 
model we use cannot give a precise answer as we have 
neglected the effects of the electrostatic self-energy of the 
free counterions. Further, as the concentration increases 
it becomes harder to distinguish the "inside" and "out- 
side" regions of the chain. We can obtain, however, a 
clear limit for the validity of these approximations, that 
also signals the onset of new phenomena. In the absence 
of multivalent counterions, the fraction of monovalent 
counterions condensed is given to a good approximation 
by 



1 



In 



2Bj\nN 



(15) 



At the onset of the dilute regime = b^^N~'^, we recover 
the Manning formula 1 f — l/Bj. Further increase on 
the concentration leads to a situation in which the frac- 
tion condensed can reach a value very close to 1, so that 
the effective charge is almost zero and cannot extend the 
chain. After this point the distinction between the col- 
lapsed and non-collapsed states is much less clearly de- 
fined. The chain cannot be extended if the electrostatic 
free energy per monomer smaller than 1/cbT (for flexible 
polymers). This leads then to the following criteria: 

-lnb^{(j> + s) {BjlnNy/^. (16) 



D. Monovalent Collapse 

As mentioned in the introduction, it is expected that 
monovalent counterions on themselves can collapse the 
chain even in dilute conditions, if the temperature is suf- 
ficiently low. Our model gives a simple criteria for the 
occurrence of this collapse. The free energy for the ex- 
tended state can be approximated as 



(lnb3(^ + g))2 
4:Bj\nN 



(17) 



while the energy for the collapsed state is simply approx- 
imated by 



-Bj. 



This results in a transition occurring at 
-ln6^(0 + s) - BjiliiN)^^^. 



(18) 



(19) 



We remark that this equation differs from Eq. (p^, by 
a Bj^/^ factor. This difference implies that the collapse 
of the chain by its natural monovalent counterions occur 
within the dilute regime. 
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V. RESULTS 



Solution of the equilibrium equations between the two 
phases produces a phase diagram that is valid at low 
concentrations of monomers. We now present results for 
valences z = 2,3, and Bj — 2,3,4,5. Other parameters 
take values of = 10^, b'^s — 10~^^, while the data cov- 
ers monomer concentrations ranging from b^cf) = 10^^*, 
to 10~*. These results are presented in Fig. 1-4. Figures 
1 and 3 present results for the concentration m, while 2 
and 4 show the required ratio to/0. 
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In both cases, z = 2, and z = 3, it is clear that as the 
density of monomer is decreased, the required amount of 
monovalent counterions required for the collapse falls in 
the range prescribed by the limiting relation Eq. ([l4|), 
and the ratio is always somewhat smaller than l/z. 

A. Limits of Applicability 

In this subsections we collect the main limitations of 
our approach. The concentration of multivalent ions at 
the extended-to-coUapse transition m* appears to be lin- 
ear with the concentration of charged monomers 4> at 
highly dilute solutions, in agreement with experimental 
results 1^. As 4> increases, however, to* decreases. Ex- 
perimental results indicate the opposite behavior [^,^. 
This disagreement has two sources: we are neglecting the 
interchain interactions and also the possibility of other 
conformations . 

In the collapsed state the chains are nearly neutral 
so one can neglect the chain-chain interactions provided 
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the concentration is still much lower than the overlap 
concentration. However, this is not the case of other 
non-collapsed conformations which have higher effective 
charge, such as the extended state. Our model takes into 
account only two states (the stretched and collapsed) of 
non-interacting chains, so that our results are valid only 
for highly dilute solution of polyelectrolytes, as has been 



roughly established in Eq. (16). In a more dense regime 



the electrostatic free energy of the free counterions be- 
comes relevant while here it has been neglected. Also, 
our results are only reliable in low ionic strength so- 
lutions of monovalent ions and at low concentration of 
multivalent ions. When the concentration of multivalent 
ions increases there will be an excess of free multivalent 
counterions so their electrostatic contribution cannot be 
neglected. We have also not taken into account the ef- 
fects of the solvent quality. In poor solvents, for exmaple, 
it is know that a "necklace" structure might be more fa- 
vorable than the stretched conformation [p7|-p^ . 

VI. CONCLUSIONS 

A collapsed flexible polyelectrolyte in low ionic 
strength solutions has a nearly zero effective charge. The 
charge of the monomers is compensated by condensed 
counterions. The main contribution to the electrostatic 
free energy of the collapsed polyelelctrolyte is propor- 
tional to the number of compensated monomer charges 
N, and the internal state of the collapsed molecule is 
similar to an amorphous ionic solid. There is a strong 
electrostatic free energy reduction of the collapsed state 
when the valence of the condensed counterions increases. 
Therefore, when multivalent counterions are added at a 
concentration m less than the concentration of charged 
monomers (j), all of the multivalent counterions are con- 
densed, replacing the monovalent condensed counterions. 
Since the multivalent ions reduce the electrostatic free 
energy, a first order transition occurs from extended to 
collapsed chain conformations. 
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